Persistent current in a mesoscopic ring with diffuse surface scattering 
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The persistent current in a clean mesoscopic ring with ballistic electron motion is calculated. The 
particle dynamics inside a ring is assumed to be chaotic due to scattering at the surface irregularities 
of atomic size. This allows one to use the so-called "ballistic" supersymmetric a model for calculation 
of the two-level correlation function in the presence of a nonzero magnetic flux. 
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Since the pioneering work of Biittiker, Imry and Lan- 
dauer |jj , persistent currents in small metallic rings have 
been a subject of great theoretical and experimental in- 
terest. While earlier observations were made on disor- 
dered samples with diffusive electron dynamics [Q, more 
recent experiments measured the persistent currents in 
high mobility semiconductor heterostructures with the 
elastic mean free path I ~ 1.3L (L is the system size) [||. 
The semiclassical electron dynamics in such clean sam- 
ples is ballistic and can be either regular (integrable) or 
chaotic. The former possibility seems to be rather excep- 
tional, since any deviation of the sample shape from a 
perfect annulus, however small, breaks down the integra- 
bility of a system. The purpose of this paper is to calcu- 
late the persistent current in a ballistic ring, in which a 
bulk disorder is absent but the electron dynamics is nev- 
ertheless chaotic due to the multiple surface scattering. 
The electron-electron interaction is neglected. 

According to Ref. j| , the average current for a canon- 
ical ensemble (i.e., for a fixed number of electrons in a 
sample) is given by the following formula: 
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(we use the units in which K = c = 1). Here s is the 
spin degeneracy, A = (poV)^ 1 is the mean level spac- 
ing in the system, p is the average density of states, 
V is the system volume, n(e) is the Fermi distribution 
function, ^(61,62;$) is the dimensionless two-level cor- 
relation function at a nonzero magnetic flux which 
depends on the energy difference e 2 — ei = 
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(Sp(E, $) = p(E, $) — po is the deviation of the one- 
particle density of states from its average value). We 
restrict our analysis to the limit of T 3> A, where the 
formula (Q) is valid. 

We assume that the sample has the shape of a pla- 
nar coaxial ring with outer and inner radii R\ and R 2 , 
respectively, threaded by a solenoid carrying a flux $ 
(see Fig. fy). We consider here only narrow rings with 
d = Ri - Rz <C R = (Ri + R 2 )/2. For this geometry, 
the number of transverse channels N = mv^djlix , the 



average density of states po — m/2ir and the mean level 
spacing A = (mRd)^ 1 . 

In order to calculate the two-level correlation function 
(^), we use the supersymmetric a model, which has be- 
come a powerful tool in the theory of disordered metals 
H and has been adapted to the description of the classi- 
cally chaotic systems as well HQ . Here it is appropriate 
to emphasize that in a clean ring, where a natural ensem- 
ble is absent but the dynamics is classically chaotic, the 
averaging in Eq. (^) must be performed over a wide en- 
ergy band , in contrast to the case of a disordered ring, 
where the angular brackets in Eq. (^) imply averaging 
over different realizations of the random potential. If the 
shape of a sample is highly symmetric and its surface is 
smooth on the atomic scale, then the specular boundary 
conditions commonly used in chaotic billiards give rise to 
integrability of the system, and the whole approach based 
on using the supersymmetric a model fails (this case was 
studied in Ref. ||, using the semiclassical trace formu- 
las). In order to make the particle dynamics chaotic, 
one can, for instance, slightly deform the shape of the 
billiard to break its perfect rotational symmetry [p|,^0[. 
However, in this approach, the level correlations can be 
calculated only numerically. Another way to achieve the 
chaotic regime, at the same time preserving the macro- 
scopic symmetry of a sample, is to assume that each act 
of the surface reflection is stochastic itself, i.e., the in- 
cident particle gets reflected in some random direction 
at the surface. This model is commonly referred to as 
the diffuse reflection and applies to the surfaces which 
are rough on the atomic scale, which seems to be quite 
reasonable physical assumption jll[] . After several reflec- 
tions at the walls, the dynamics of an electron becomes 
fully chaotic. 

At lo ^> A, the two-level correlation function (||) can be 
calculated perturbatively, using the "ballistic" version of 
the supersymmetric a model generalized to the presence 
of a nonzero magnetic flux |L2| : 
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where Ai($) are the eigenvalues of the (Cooperon) Liou- 
ville operator 
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inside the Aharonov-Bohm billiard (n is the direction of 
momentum). The region of small frequencies to <C A, 
where the perturbative approach fails and the level cor- 
relations are described by the universal formulas of the 
random matrix theory 13 1, lies beyond the limits of ap- 
plicability of the thermodynamic approach to the descrip- 
tion of persistent currents. The expression ([}]) is a direct 
analog of the Altshulcr-Shklovskii spectral function for 
diffusive systems E3. Rewriting the Fermi distribution 
functions in Eq. (II) as Matsubara sums and integrating 
over energies, we end up with the following expression: 
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where uj n = 2ixnT . The sum over n on the right-hand 
side is convergent due to the presence of the differentia- 
tion over flux. 

The spectrum of the Liouville operator is determined 
by the eigenvalue equation 
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with some boundary conditions at the surfaces of the 
ring. Due to the similarity of Eq. (^) to the Boltz- 
mann kinetic equation, /(r, n) having the meaning of 
the classical distribution function, the boundary con- 
ditions at the diffusely reflecting surface r = R with 
the outward normal N can be imposed by analogy with 
the classical kinetic theory. The distribution function 
of reflected particles can be represented as /(R, n) = 
p/o(R) + (l-p)/(R, n) @, where < p < 1 is "the dif- 
fuseness coefficient" , f is an isotropic distribution func- 
tion, and ri = n 2(nN)N is the direction of specular 
reflection. In this paper, we consider an isotropic dif- 
fuse scattering with p = 1, corresponding to the limit 
of "strong chaos" (for the discussion of applicability of 
this model to real experimental samples, see below) and 
the boundary condition, which follows from the particle 
number conservation, takes the form 

-/(Ri,n)| (nNi)<0 = / dn'tn'NO/^.n'). (7) 

n J(n'Ni)>0 

Here J dn = J dcj)/2ir ((f> is the angle between the direc- 
tion of momentum of incident particles and the outward 
normal N), and i — 1,2 correspond to the outer and 
inner surfaces of the ring. The distribution function of 
reflected particles on the left-hand side of Eq. (fj]) does 
not depend on n. A similar approach was used in Ref. 
jl6| for calculation of the corrections to the universal level 
correlations in a two-dimensional disk without magnetic 
flux. 

Since the magnetic field is absent inside the ring, the 
trajectory of a particle between collisions with the walls 
is a straight line. Equation (|TJ) can be solved along the 
trajectory 
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This expression allows one to establish a rela- 
tion between /(Ri, n)| (nNl)>0 and /(Ri, n)| (nNl)<0 or 
/(R-2,n)| (nN2 ) <0 , and also between /(R 2 , n)| (nN2)>0 and 
/(Ri, n)|( nNl ) <0 in Eq. (]?]) and obtain a rather cumber- 
some algebraic equation for the eigenvalues of the Liou- 
ville operator in a ring of arbitrary width. Fortunately, 
in the case of a narrow ring, the problem can be consider- 
ably simplified, since we can replace our annular billiard 
by a strip of length L = 2nR and width d such that 
8 = d/L <C 1 (see Fig. ||). In addition, the vector poten- 
tial can be put constant inside the sample: Ag = <1>/L 
(8 is the polar angle in real space). This simplifica- 
tion implies that the contribution from the trajectories 
connecting two points at the outer wall is neglected (it 
can be checked that this contribution is indeed small at 
i5< 1). However, there is an important property of the 
annular geometry which should be taken into account, 
namely the finiteness of the flight length i(</>) between 
successive collisions with the walls (I «C L). This feature 
can be restored in the strip billiard if to assume that 
there exists the maximum scattering angle (f>o such that 
sin</>o = R2/R1 ~ 1 — 2tt5. 

Let x = R8, then we obtain, from Eq. (0) and Fig. 0: 



h,>{x,4>) = h,<i x ~ dtan0) 

xex p(^^) ex p( 
h,>(x,4>) = fx,<{x + dt&ncj)) 



9 ■ 2$ dtan0 



(9) 



xex p(^^) ex p(- 



27ri 



2$ d tan <; 
¥0 T~ 



Here /i,>(<)(x, (j>) = /(R i5 n)| (nNi)>0(<0 ) and <J> = 
2irhc/e is the flux quantum. It is convenient to ex- 
pand the functions ft in the Fourier series: fi.<(x) = 
J2 q fi,q e%qx ■ Since f(x + L, n) = f(x, n), the wave num- 
ber is quantized: q = m(2Tr/L), where m = 0,±1,.... 
Substitution of the expressions (^) in Eq. (^) results 
in the following equation for the eigenvalues A m .fc( < I > ) = 
(vp I 'd)z m ^k(<&) of the Liouville operator (|J): 



F m (z,u) = A 2 m (z,u) -1 = 0. 
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Here v = 2$/$ and 
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To guarantee correct normalization of the distribution 
function of reflected particles, the prefactor sin -1 <j>o has 
been included in the definition of A m . 

The eigenvalues of the Liouville operator at fixed m 
are labeled by k — 0,1,2,..., and have complex conju- 
gated partners for all (m, fc). At $ = 0, Eq. (10|) has the 
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solution zo.o — corresponding to the equilibrium distri- 
bution function, which does not depend on n and r. It 
is clear from Eq. ( |Tl| ) that the eigenvalues \ m ,k(&) have 
the following property: \ m , k (® -I- $o/2) = A m+ i jfc ($) so 
that any quantity which can be represented in the form 
Em fc ' l (Am,fc( < I > )), where h(X) is some function, must be 
periodic function of $ with the period of half the flux 
quantum $0. Note that, since the eigenvalues of the Li- 
ouville operator physically correspond to the relaxation 
rates of different harmonics of a nonequilibrium classi- 
cal distribution function, the real parts of all \ m> k are 
positive. 

Using the representation of the sum over i — (to, k) in 
Eqs. (|) and ([§ as an integral over a contour C enclosing 
all zeros of the function F m [v, z) in the complex plane of 
z: 
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where Jo = ev F / L is the current carried by a single elec- 
tron state in an ideal one-dimensional ring, and F m {z, v) 
is given by Eq. (|lC)|). 

Due to the existence of different energy scales in the 
system, the temperature dependence of the persistent 
current is characterized by several distinct regimes. The 
smallest energy scale is the mean level spacing A, which 
also limits the applicability of the thermodynamic ap- 
proach itself. Two other scales are given by the in- 
verse times i^ 1 = Vf/L — NA and t^ 1 — vp/d — 
NA/6. It follows from Eqs. @ and (|ll|) that at 
T 3> t^ 1 the persistent current is exponentially small: 
7~ I exp(-T//VA). 

In a multichannel ring, there also exists yet another 
energy scale A <C E c <C NA, whose origin can be most 
easily understood if to return to Eq. (||) . Using the iden- 
tity a~ 2 = J Q dye~ ay y and calculating the sum over n, 
we end up with the following expression 
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At sufficiently low temperatures, the dominant contribu- 
tion comes from Ao.o( < I > ), which is the solution of Eq. 
( |To| ) with smallest real part. At $ — > 0, the energy 
E c coincides with Ao.i( < I > = 0), and at T <C E c the 
sum over (m, fc) in Eq. ( |l4| ) can be replaced by its 
value at m = k = 0. It follows from Eq. ( |lo| ) that 
zofiiv) = 2tt8 2 \ \t\5\v 2 at v —y 0, and we obtain: 
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Thus the persistent current and also the orbital mag- 
netic moment M = ttR 2 I of a ballistic ring exhibit a 
Curie- type response on magnetic flux at A <C T <C E c . 

The energy E c is associated with the relaxation time 
t c = E~ l of any initially nonequilibrium state of the 
system to a spatially uniform and isotropic distribution 
(the ergodic time) . In the systems with diffusive electron 
dynamics, t c = L 2 j D (L is the system size, D is the dif- 
fusion coefficient) and E c is called Thouless energy. In 
general, it is not so easy, however, to say what the er- 
godic time is in ballistic systems 0. In our case, it can 
be estimated from simple physical considerations, with- 
out solving Eq. ([n]) explicitly, as follows. The typical 
flight time tf between two reflections (see Fig. 0) is 
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Since the particle completely "forgets" the initial di- 
rection of its momentum after each collision with the 
walls, the characteristic time for establishing an isotropic 
momentum distribution is of the order of tf. In con- 
trast, the characteristic time of filling uniformly all avail- 
able configuration space is much longer. Indeed, the 
typical displacement between two collisions is given by 
\J {(5x) 2 ) ~ Vptf ~ L<5 3 / 4 . According to the central 
limit theorem, the probability to find a particle at the 
distance Ax = Yli=i $ x i after M ^> 1 collisions obeys 
th e Gauss ian d istribution with the standard deviation 
^/((Ax) 2 ) = ^/M{{8x) 2 ). Substituting here M ~ t/t f 
and using Eq. (|l6|), we obtain 
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Thus the motion of an electron along the circumference 
of the ring is in fact diffusive at t 3> tf, and 
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In a narrow multichannel ring the hierarchy of the charac- 
teristic energy scales looks as follows: A <C E c <C t^ 1 <C 
tj 1 <C i^ 1 . Note that all our results are valid if one can 
neglect the deflection of the electron trajectories inside 
the ring by an external magnetic field, which is the case if 
the inverse cyclotron radius R~ l — (eB/mvF) is smaller 
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than the inverse typical flight length (vptf) 1 . This con- 
dition, rewritten as <&/<&o N5~ 7 / 4 , is always satisfied 
in narrow rings. 

In order to facilitate comparison of our results with 
the experimental data, let us rewrite Eq. fll5| ) in a differ- 
ent form, using the identity A/T = N~ 1 (Lt / L), where 
Lt = vf/T is the length scale associated with temper- 
ature. In the experiment of Ref. ||, 5 ~ 0.1, N ~ 4, 
Lt/L ~ 5, so that T^A. Due to the presence of the 
factor 5 In 5 <C 1 in Eq. ([[5]) , the predicted current turns 
out to be smaller than the experimentally observed (and 
also than the theoretically calculated for the case of spec- 
ular reflection ||). This discrepancy can be attributed 
to the fact that, because of the low density of carriers, 
the Fermi wavelength greatly exceeds the size of the sur- 
face irregularities, so that a considerable fraction of par- 
ticles gets reflected specularly rather than diffusely (i.e., 
p < 1), and one should describe the semiclassical dynam- 
ics in the experimental conditions of Ref. || as "weakly 
chaotic" . 

To summarize, we calculated the persistent current in 
a small clean metal ring, in which the electron dynam- 
ics is chaotic due to the stochastic surface scattering. A 
general analytical expression for the persistent current 
is derived in the limit of "strong chaos", and a Curie- 
type orbital magnetic response on a small external flux 
is predicted at A < T < E c . 
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FIG. 2. The transformation of a ring into a narrow strip 
of width d = Ri - R 2 and length L = 2-kR = 7t(Ri + R 2 ). 
A — &/2nR is the vector potential created by the solenoid. 
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